Large temperature dependence of the number of carriers in Co-doped BaFe 2 As 2 
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Using angle-resolved photoemission spectroscopy, we study the evolution of the number of carriers 
in Ba(Fei_ x Co x )2As2 as a function of Co content and temperature. We show that there is a In- 
dependent energy shift compared to density functional calculations, which is large at low Co contents 
and low temperatures and reduces the volume of hole and electron pockets by a factor 2. This k-shift 
becomes negligible at high Co content and could be due to interband charge or spin fluctuations. We 
further reveal that the bands shift with temperature, changing significantly the number of carriers 
they contain (up to 50%). We explain this evolution by thermal excitations of carriers among the 
narrow bands, possibly combined with a temperature evolution of the k-dependent fluctuations. 



PACS numbers: 79.60.-i, 71.18.-y, 71.30.-h 

Since the discovery of iron based superconductors, an- 
gle resolved photoelectron spectroscopy (ARPES) has 
given valuable information about their electronic struc- 
ture [UH]. In most cases, the measured spectra are in 
broad agreement with Density Functional Theory (DFT) 
calculations, after renormalization by a factor 2 to 3. 
This range of values is in good agreement with mass en- 
hancements predicted by Dynamical Mean Field Theory 
(DMFT) [3l|4]. More unusual is a shrinking of the hole 
and electron pocket sizes compared to DFT calculations, 
observed both by ARPES and de Haas- Van Alphen ex- 
periments [5 j. In a Fermi liquid (FL) picture, this cor- 
responds to a down shift of the hole bands compensated 
by an up shift of the electron ones, as sketched in Fig. 
1(a). As it depends on k, we refer to this effect as k- 
shift in the following. Its origin is so far unclear, but 
it seems associated with the strength of correlations. In 
P-substituted BaFe2As2, dHvA measurements reported 
that the maximum k-shift corresponds to the largest ef- 
fective masses [6]. In Ru-substituted BaFe2As2, both the 
k-shift and the renormalization weaken at large doping 
content [TJ [8] . Ortenzi et at. showed that interband in- 
teractions mediated by low energy bosons, such as spin 
fluctuations, can lead to such a k-shift [9]. If this is indeed 
the case in pnictides, it contains important information 
about the interactions in these compounds and deserves 
a detailed study. 

We show here with ARPES in Co-doped BaFe2As2 
that the k-shift is suppressed for high Co content, when 
the hole pockets fill up. We reveal in addition a very 
strong temperature (T) dependence of the electronic 
structure. Up to 10% Co, the electron pockets are ex- 
panding with increasing T by an amount as large as 50%. 



A similar T dependence was reported recently for Ru- 
doped BaFe 2 As 2 [10 . At higher Co content (30%), the 
electron pockets on the contrary shrink at high T. We 
show that this can be largely explained by thermal exci- 
tations in these narrow multiband systems. We present 
both simple models and calculations using the density of 
states (DOS) obtained by DFT within the Local Density 
Approximation (LDA) and by DMFT in order to sim- 
ulate the effect of temperature within a FL framework. 
These simulations allow to account for about half the ex- 
perimental shift. Part of the expansion of the electron 
pockets could also be due to the weakening of the k-shift, 
as proposed in [11] for the model of interband interac- 
tions, or to an evolution of correlation effects. 

Single crystals of Ba(Fei_ :E Co :E )2As2 were grown using 
a FeAs self-flux method and studied in detail by trans- 
port measurements [12] . ARPES experiments were car- 
ried out at the CASSIOPEE beamline at the SOLEIL 
synchrotron, with a Scienta R4000 analyzer, an angular 
resolution of 0.2° and an energy resolution of 15 meV. 
All measurements were perfomed at 34 eV, except other- 
wise mentioned. DFT band structure calculations were 
performed using the Wien2k package [13] at the experi- 
mental structure of BaFe2As2- To model Co doping, we 
used the virtual crystal approximation (VCA), where Co 
induces a global rigid shift of the entire band structure. 
For the LDA+DMFT calculation, we used the technol- 
ogy of ref [4 , which fully takes into account dynamical 
screening effects from first principles. 

Recently, we have shown that the electron pocket 
around X is an ellipse formed by a "deep" electron band 
of d xy symmetry along the long axis and a "shallow" elec- 
tron band of d xz /d yz symmetry along the short axis [15] • 
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Figure 1. (a) Calculated LDA band structure along TX at k z =0 and x=0. Colors indicate the main orbital character. Arrows 
indicate the direction of the k-shift observed at low T in experiments, (b) Jgf measured as a function of Co content for the two 
electron bands, deep and shallow. The thin lines are guides for the eye. Solid lines are Uf obtained in VCA DFT calculations 
for k z corresponding to experimental conditions [14]. (c) Shifts compared to LDA calculations at x=0 needed to obtain the 
experimental /cf- Solid line indicates the expectations from the VCA calculations. 



In Fig. 1(b), we report kp for these two bands as a 
function of Co doping at 20 K (150 K for BaFe2As2 to 
avoid the magnetic state). In the supplementary infor- 
mation [14], we show the full dispersions used to extract 
these values. The pocket expands with Co content, as 
expected from the effective electron doping of the bands 
[T6HT9] . For small Co contents, the values are signifi- 
cantly smaller in the experiment than in the DFT VCA 
calculations, shown as solid lines. This reflects the up 
shift of the electron bands sketched in Fig. 1(a) and de- 
scribed in the introduction. Interestingly, the shift also 
increases significantly the anisotropy of the pocket near 
x=0 doping. 

In Fig. 1(c), we report as a function of Co content the 
shifts with respect to the LDA bands calculated at x=0 
needed to obtain the measured kp. As kp is uniquely 
determined by the number of carriers in one pocket, this 
shift value is independent of the band renormalization 
and well defined. We found similar shift values for dif- 
ferent k z [15]. Compared to the rigid shift found in the 
DFT VCA calculation (orange solid line), the bare shift is 
nearly 100 meV at low dopings, but it vanishes at higher 
dopings. In Table I, we quote the shift at x=0, as well as 
the renormalization extracted from the dispersion [T4] . 

It is more difficult to perform a similar analysis for 
the hole pockets, because the shifts and renormalization 
values are strongly /^-dependent (see Table I and [T4]). 
Such a deformation of the band structure is not surpris- 
ing considering that the bands must shift in opposite di- 
rections at T and X. Fig. 1(a) recalls for example that 
the d yz hole band must be degenerate at X with the d xz 
electron band. As this point is fixed in energy along k z 
by symmetry, while the hole bands shift quite strongly, 
this implies some modifications of the hole pockets k z 
dependence. We report in Fig. 1(c) the shift of the cal- 
culated d xz /d yz bands needed to describe the top of the 
hole band at k z =0. Like for electron pockets, the k-shift 
disappears quite precisely when the hole pockets become 
completely filled, near x=0.2 [16], suppressing the possi- 
bility of adjusting the number of carriers by relative shifts 



of the bands. 

The temperature dependence should give interesting 
and complementary information to explore the nature of 
the k-shift. If related to spin fluctuations, one could ex- 
pect interband interactions to weaken at high T, as the 
strength of spin fluctuations decreases [IT] . In Fig. 2, we 
show T dependent measurements from 20 to 300 K for 
8% and 30% Co contents. The full dispersion images are 
shown in the supplementary information [14 . In each 
case, we took care to check that the effect was reversible, 
as there is often some evolution with temperature cycles 
(broadening and/or loss of intensity). In Fig. 2(a) and 
2(d), we show the spectral intensity of the shallow band 
integrated in a 10 meV window around the chemical po- 
tential \i. The distance at a given temperature between 
the two maxima defines 2kp. White dashed lines indicate 
its values at low T to emphasize the evolution at higher 
T. As reported by red squares in Fig. 2(b) and 2(e), kp 
seems roughly constant up to 150K and then increases at 
8% Co and decreases at 30% Co. For the other bands, 
there is also an evolution of kp with T, which is reported 
in the same figures. For Co 30%, the hole pockets are 
filled and small electron pockets appear at T [THUS]. The 
opposite evolution of kp for the shallow electron pocket 
for the two dopings evidences that the evolution depends 
on the nature of the carriers at \i. In fact, these changes 
of kp correspond to a shift of the whole bands [14]. In 
Fig. 2(c) and 2(f), we report these shifts as a function 
of T. For Co 8%, they are remarkably similar, indicating 
that the evolution in temperature is mainly driven by a 
shift of the \i within the electronic structure. 



Table I. Shifts and effective masses with respect to LDA 
bands calculated at x=0 [14] . 
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Figure 2. Temperature dependence for 8% Co (a-c) and 30% 
Co (d-f). (a,d) Spectral weight integrated over lOmeV around 
\i for the shallow electron band. White dashed line indicate \zf 
at T=0. The last row (20 K) was measured after cooling back 
from 300K. (b,e) fcj? as a function of T for the two electron 
bands at X (shallow and deep) and for one band at Y (holes 
for x=0.08, electrons for x=0.3). (c,f) Shift as a function of 
T for the different bands. 



In Fig. 3(a), we give a simple sketch of the Co 8% 
band structure with realistic parameters : a hole band 
whose top is 10 meV above Ep (in experiment, it is even 
below Ep for a large range of k 2 [14j) and an electron 
band whose bottom is 50 meV below Ep, as for the shal- 
low electron band. We assume a constant DOS p(E) and 
choose it so that both bands contain the same number 
of carriers no at T=0. At 300K, the Fermi edge (thin 
line) extends over the top or bottom of the bands. If the 
chemical potential p did not move, the number of carriers 
n = J p(E)f(E) dE would be quite different from T=0, 
namely 1.07no for electrons and 2.33no for holes in this 
example. To keep the number of carriers constant, p has 
to shift to a higher position that will reduce the number 
of excited carriers. As the imbalance for holes is much 
larger than for electrons, p stabilizes a different value of 
n compared to T=0. A simple calculation shows that the 
new number of carriers would be 1.38no [Fig. 3(b)] and 
p should shift by 17 meV at 300 K to ensure n e =n/ l [blue 
line in Fig. 3(d)]. Note that, despite the fact that the 
hole bands are now completely below Ep, as in our mea- 
surement, they do contain carriers. Similarly, the fact 
that kp at Co 8% decreases for holes and increases for 




Figure 3. (a) Thick lines : DOS considered for the 2 bands 
model. The Fermi edge is sketched at 300K, dashed areas indi- 
cate the occupied parts for holes and electrons before shifting 
the chemical potential fi. Arrows indicate the directions in 
which the bands will shift with increasing T with respect to 
/i. (b) Number of carriers (holes or electrons) obtained as a 
function of temperature for the 2 bands model, (c) DOS cal- 
culated in LDA, LDA+DMFT and LDA+DMFT including a 
k-shift. The Fermi levels for x=0 and x=0.3 are indicated by 
vertical lines, (d) Shift of \i obtained within band theory from 
the DOS of various models (see text and caption). 



electrons does not violate the conservation of the num- 
ber of carriers. Although very rough, this model gives the 
right order of magnitude for the observed shift, support- 
ing the idea that this effect plays a leading role. To go 
further, we now consider in Fig. 3(c) a realistic DOS for 
BaFe2As2- As in our simple model, it is the occurrence of 
fine structures within 4/c#T of the Fermi level that will 
induce a shift of p. Indeed, within FL theory, one ex- 
pects bp{E) oc T 2 p'(E)/p(E) [20 . In the LDA DOS, the 
peaks are too far from Ep to induce any sizable shift and 
p remains constant in Fig. 3(d). As a next step, we use 
the DMFT spectral function calculated at 150 K as an 
approximation for the DOS. The renormalization of the 
bands bring the peaks closer to Ep and a sizable shift is 
detectable at 300K in Fig. 3(d). This shift even increases 
when the bands extrema are moved closer to Ep by in- 
troducing arbitrarily a k-dependent shift simulating that 
observed experimentally [21]. For x=0.3, Ep moves to 
a region of opposite slopes inducing a change of sign in 
the shift, as observed in Fig. 2(f). Even if the values are 
somewhat smaller than in the experiment, this captures 
the qualitative behavior very well. As will be discussed 
in the conclusions, a T dependence of the k-shift could 
also contribute to these effects. 

In Fig. 4, we finally compute the number of electrons 
as a function of temperature and doping. Thanks to the 
knowledge of the 3D shape of the electron pockets [15], 
this can be reliably obtained by combining the kp mea- 
sured for the deep and shallow bands and integrating 
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over k z [22]. In the right part of the figure, we report 
the change of the number of electrons as a function of 
T. For 8% Co, this leads to an increase by 50% between 
low and high temperature. This order of magnitude is 
in good agreement with the simple model of Fig. 3. For 
BaFe2As2 above the magnetic transition, the effect is also 
present and even slightly larger. This rules out that the 
anomalous temperature dependence could be due to per- 
turbations associated with the presence of Co. For 15% 
Co, the effect is somewhat reduced. For 30% Co, the 
number of electrons at X decreases with increasing T [see 
kp in Fig. 2(e)], but it is difficult to evaluate the num- 
ber of electrons at T, so that we only indicate the low T 
value [23 . The left part of Fig. 4 displays the number of 
electrons as a function of Co content. This emphasizes 
the amplitude of the change in temperature compared to 
that with doping. It also shows that the T dependence 
tends to reduce the deficit of electrons present at low T 
compared to the calculation. For 30% Co, the hole pock- 
ets are completely filled and the number of electrons is 
almost that expected by stoichiometry. 

Interestingly, this T dependence of the number of elec- 
trons is strikingly similar to the one deduced from Hall 
measurements [12] , although a bit smaller quantitatively. 
This strongly supports the interpretation of both ARPES 
and transport based on a change of the number of elec- 
trons with T. As a result, much care should be taken to 
extract the genuine T dependence of the scattering times 
from resistivity measurements. Another quantity that 
exhibits a quite similar and anomalous T dependence is 
the spin susceptibility as measured by NMR for exam- 
ple [24]. In this case however, \ increases with temper- 
ature but decreases with electron doping. Of course, x 
should not solely depend on the number of electrons, but 
also on the number of holes, for example. In any case, it 
will be important to consider this T dependent number 
of carriers to reinvestigate these properties [25] . 

Finally, our study reveals peculiar and so far largely 
unexplored consequences of the multiband nature of iron- 
based superconductors. First, there is a shift compared 
to band calculations that is k-dependent and leads to the 
shrinking of the Fermi Surface for low Co content. This k- 
shift has a dramatic effect as it reduces the total number 
of carriers by a factor of 2 for x=0. It also automatically 
leads to more elliptic electron pockets and more 3D hole 
pockets than in the LDA band calculation, as we observe 
by ARPES. These k-shifts weaken with doping, as the 
hole bands fill up, which is consistent with the picture 
based on interband couplings [9j [11] . It would be very 
interesting to obtain a full theoretical description of this 
effect. It has not been seen in DMFT based electronic 
structure calculations up to now, possibly because it re- 
quires an extremely high energy precision to be resolved 
(of the order of 10 meV). Alternatively, one may spec- 
ulate that it is due to non-local self-energy effects not 
included at the DMFT level. A very recent GW study 
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Figure 4. Number of electrons as a function of Co content for 
different temperatures (left part) and as a function of T for 
different Co content (right part). The lowest T is shown as 
filled symbol. The orange thick line indicates the number of 
electrons expected from LDA calculation. 

may be interpreted in this sense [26]. Independently of 
their origin, an important consequence of these k-shifts 
is that they make the whole band structure much more 
sensitive to temperature than what would be expected in 
most calculations. Indeed the most striking finding of our 
work is a strong dependence of the electronic structure on 
temperature, confirming and generalizing the finding of 
a previous ARPES study [ID] . The resulting large T de- 
pendence of the number of carriers must impact all other 
electronic properties. There is already good evidence that 
transport probes a similar evolution [12]. We show that 
thermal excitations within the FL picture partly explain 
the effect. A full theoretical description will probably 
requires a complete understanding of the origin of the 
k-shift, so that its T dependence can be controlled. 
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